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ABSTRACT 

We use Monte Carlo techniques to simulate the ability of future large high-redshift galaxy 
surveys to measure the temporal evolution of the dark energy equation-of-state w{z), using the 
baryonic acoustic oscillations in the clustering power spectrum as a 'standard ruler'. Our analysis 
only utilizes the oscillatory component of the power spectrum and not its overall shape, which 
is potentially susceptible to broadband tilts induced by a host of model-dependent systematic 
effects. Our results are therefore robust and conservative. We show that baryon oscillation 
constraints can be thought of, to high accuracy, as a direct probe of the distance — redshift and 
expansion rate — redshift relations where distances are measured in units of the sound horizon. 
Distance precisions of 1% are obtainable for a fiducial redshift survey covering 10,000 deg 2 and 
redshift range 0.5 < z < 3.5. If the dark energy is further characterized by w(z) = wo + W\Z 
(with a cut-off in the evolving term at z = 2), we can then measure the parameters wq and w\ 
with a precision exceeding current knowledge by a factor of ten: ler accuracies Au>o rj 0.03 and 
Aw\ 0.06 are obtainable (assuming a flat universe and that the other cosmological parameters 
fi m and h could be measured independently to a precision of ±0.01 by combinations of future 
CMB and other experiments). We quantify how this performance degrades with redshift/areal 
coverage and knowledge of f2 m and h, and discuss realistic observational prospects for such large- 
scale spectroscopic redshift surveys, with a variety of diverse techniques. We also quantify how 
large photometric redshift imaging surveys could be utilized to produce measurements of (wq, wi) 
with the baryonic oscillation method which may be competitive in the short term. 

Subject headings: cosmological parameters — large-scale structure of universe — surveys 



1. Introduction 

In the current standard cosmological model, 
baryonic matter and cold dark matter together 
contribute only about a third of the total energy 
density of the Universe. One of the most impor- 
tant puzzles in cosmology is to account for the 
remaining two-thirds of the energy, which is re- 
quired to render the Universe approximately spa- 



tially flat, as demanded by recent observations 
of the Cosmic Microwave Background (CMB; de 
Bernardis et al. 2000). 

The existing set of cosmological data - con- 
sisting principally of measurements of the CMB, 
of the local clustering of galaxies, and of distant 
supernovae - can be understood by invoking the 
existence of the 'cosmological constant' A orig- 
inally envisaged by Einstein, such that it con- 
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tributes a present-day fractional energy density 
n A = 0.73 ± 0.04 (Spergel et al. 2003). The re- 
markable consequence of this model is that A acts 
as a 'repulsive gravity', driving the rate of cosmic 
expansion into a phase of acceleration. Equally 
surprisingly this acceleration has been observed 
reasonably directly by an apparent dimming of 
distant supernovae (Riess et al. 1998; Perlmutter 
et al. 1999). 

The cosmological constant component is natu- 
rally attributed to the inherent energy density of 
the vacuum. However, the 'expected' quantum- 
mechanical Planck energy density is larger than 
that required to account for the accelerating rate 
of cosmic expansion by an exceptionally large 
dimcnsionless factor of order c 5 G~ 1 H~ 1 Hq 2 <~ 
10 122 . This profound difficulty has motivated the 
development of alternative models for the 'dark 
energy' - i.e. the causative agent of accelerat- 
ing expansion. Many of these models, such as 
'quintessence' (Ratra & Peebles 1988), feature a 
dynamic component of dark energy whose prop- 
erties evolve with time (e.g. a rolling scalar field 
4>). These predictions are commonly character- 
ized in terms of the dark energy equation-of-state 
w{z) = P/p, relating its pressure P to its energy 
density p (in units where the speed-of- light c = 1). 
For the cosmological constant model, w{z) = — 1 
at all epochs. 

These competing models for the dark energy 
are essentially untested, because the current cos- 
mological dataset is not adequate for delineat- 
ing variations in the function w(z) with redshift, 
which is the essential requirement for distinguish- 
ing quintessential cosmologies from a cosmolog- 
ical constant. New experiments are demanded, 
which must be able to recover the characteristics 
of dark energy with unprecedented precision, in- 
cluding any cosmic evolution. The study of the 
nature of dark energy is the current frontier of ob- 
servational cosmology, and has been widely iden- 
tified as having profound importance for physics 
as a whole (Turner et al. 2003). 

Given the fundamental importance of accel- 
erating cosmic expansion and the possibility of 
confounding systematic effects in the supernova 
data, other precision probes of the dark energy 
model are clearly desirable. In Blake & Glaze- 
brook (2003), hereafter Paper I, we suggested that 
the small-amplitude 'baryonic oscillations', which 



should be present in the power spectrum of the 
galaxy distribution on large scales (> 30 Mpc), 
could be used as a 'standard cosmological ruler' to 
measure the properties of dark energy as a func- 
tion of cosmic epoch, provided that a sufficiently- 
large high- redshift [z > 1) galaxy survey was avail- 
able (see also Seo & Eisenstein 2003; Linder 2003; 
Hu & Haiman 2003). The characteristic sinusoidal 
ruler scale encoded by the baryonic oscillations is 
the sound horizon at recombination, denoted s. 
This length scale is set by straight-forward linear 
physics in the early Universe and its value is deter- 
mined principally by the physical matter density 
(p m cx VL m h 2 ). The cosmological uncertainty in 
this parameter combination is ameliorated by an 
advantageous cancellation of pre-factors of ^l m h 2 
in the ratio of low-rcdshift distances to s (see 
Paper I). The residual dependence of s on O m 
and other cosmological parameters is small (Eisen- 
stein & White 2004) , rendering observations of the 
acoustic oscillations a powerful geometric probe of 
the cosmological model. 

This acoustic signature has recently been iden- 
tified at low redshift in the distribution of Lumi- 
nous Red Galaxies in the Sloan Digital Sky Sur- 
vey (Eisenstein et al. 2005) (see also Cole et al. 
2005). Although these data are insufficient for pre- 
cise measurements of the dark energy, this analysis 
represents a striking validation of the technique. 
The challenge now is to create larger and deeper 
surveys. In Paper I we demonstrated that, given 
a galaxy redshift survey at z <~ 1 mapping a total 
cosmic volume several times greater than that of 
the Sloan main spectroscopic survey in the local 
Universe (we define Vsioan = 2 x 10 8 h~ 3 Mpc 3 ), 
the equation-of-state of the dark energy could be 
recovered to a precision Aw 0.1 (assuming a 
model in which w is constant). The precision of 
this experiment scales with cosmic volume V in 
a predictable manner (roughly in accordance with 
l/W) and it is not unfeasible to imagine an ul- 
timate 'all-sky' high-redshift spectroscopic survey 
within ~ 20 years. 

We believe that this standard ruler technique 
would powerfully complement proposed future su- 
pernova searches such as the SNAP project (e.g. 
Aldcring et al. 2002, http : //snap . lbl . gov) per- 
mitting, for example, direct tests of the 'reci- 
procity relation' which predicts that the true lu- 
minosity distance cIl{z) (measured by a standard 
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candle) is exactly the same as the angular di- 
ameter distance <1a{z) (measured by a standard 
ruler), up to a factor of (1 + z) 2 (Bassett & Kunz 
2004). Furthermore, we argue that it is not yet 
conclusively proven that the dimming of super- 
novae caused by cosmic acceleration can be distin- 
guished with sufficient accuracy from other possi- 
ble systematic effects, such as changes in the in- 
trinsic properties of supernovae with galactic en- 
vironment (e.g. metallicity) , dust extinction ef- 
fects, population drift and the difficulties of sub- 
percent-level photometric calibration (including 
K-corrections) across wide wavelength ranges. It 
is therefore important to pursue alternative precise 
high-rcdshift probes of the cosmological model. A 
particular advantage of the baryonic oscillations 
method is that it is probably free of major sys- 
tematic errors, assuming that galaxy biasing on 
large scales is not pathological. 

In Paper I we developed a Monte Carlo, semi- 
empirical approach to transforming synthetic 
galaxy rcdshift surveys into constraints on the 
value of w. This contrasts with, and comple- 
ments, other analytical approaches to the problem 
such as those using Fisher-Information techniques 
(e.g. Seo & Eiscnstcin 2003; Hu & Haiman 2003). 
In this present study we extend and refine our 
methodology to the more general case where we 
allow the equation-of-state of dark energy to have 
a dependence on redshift, w{z) = w + W\z, and 
we recover joint constraints upon (u> ,wi). Fur- 
thermore, we apply the standard ruler indepen- 
dently to the radial and tangential components of 
the power spectrum, which has the effect of pro- 
ducing separate measurements of the co-moving 
angular-diameter distance to the effective redshift 
of each survey slice, x(z), and the rate of change 
of this quantity with redshift, x'{z) = dx/dz. We 
note for clarity that in a flat Universe: 

x(z) =D A (z){l + z) 

x\z) = c/H{z) 

where Da is the physical angular diameter dis- 
tance and H(z) is the Hubble factor (Universal 
expansion rate) at redshift z. We also examine in 
more detail the effect of uncertainties in cosmo- 
logical parameters such as the matter density and 
the Hubble parameter. In addition, we present an 
in-depth discussion of the observational require- 
ments and prospects for realistic surveys, based 



upon both spectroscopic redshifts and photomet- 
ric redshifts. 

The plan of this paper is as follows: in Sec- 
tion 2 we give a very detailed description of our 
methodology, greatly expanding on Paper I, and 
the approximations we made to allow us to sim- 
ulate a range of large surveys, and in Section 3 
we quantify the size of redshift survey required 
to detect the oscillatory component of the power 
spectrum. The recovered constraints on (w ,wi) 
for various simulated surveys are presented in Sec- 
tion 4 for spectroscopic surveys and in Section 5 
for photometric- redshift surveys. Finally, in an 
Appendix we present the results of a very large 
computation designed to test the effects of our ap- 
proximations. 

2. Monte Carlo Methodology 

Our methodology for simulating future galaxy 
redshift surveys and assessing their efficacy for 
measuring acoustic oscillations was summarized in 
Paper I. In this Section we provide a more detailed 
account of our procedures. In addition, we have 
implemented various extensions to the methodol- 
ogy of Paper I: 

1. We now fit separate acoustic oscillation 
scales in the tangential and radial directions. 
In Paper I we fitted to the angle-averaged 
power spectrum, effectively assuming that 
the shift in the apparent radial and trans- 
verse scales as the cosmology was perturbed 
about the fiducial value were the same. This 
is only approximately true for z ~ 1 and 
breaks down at low and high rcdshift (see 
Figure 5 of Paper I). In the new approach we 
make use of the different rcdshift dependen- 
cies of radial and transverse scales to pro- 
vide extra cosmological constraints increas- 
ing signaknoise. Specifically, the tangential 
scale is controlled by the co-moving distance 
to the effective redshift of the survey, x(z), 
and the radial scale is determined by the 
rate of change of this quantity with redshift, 
x'(z) = dx/dz = c/H(z), where H(z) is the 
Hubble constant measured by an observer 
at redshift z. This is useful because H(z) is 
directly sensitive to the dark energy density 
at the redshift in question. 
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2. We allow the equation-of-state of dark en- 
ergy w(z) to have a dependence on redshift, 
w(z) = wo + wiz, and we recover joint con- 
straints upon (wq,wi). We do not claim 
that this expression faithfully describes dark 
energy in the real Universe. In particular, 
models with wi > become unphysical at 
high redshift unless we impose a cut-off for 
the evolving term: we assume that w(z > 
2 C ut) = w a + wiz cut where z cut = 2, and we 
ensure that matter dominates at high red- 
shift, i.e. wi < {—Wo)/z cu t. However, usage 
of the equation w(z) = wq + ui\z facilitates 
comparison with other work such as Seo 
& Eisenstein (2003) who use the same pa- 
rameterization, empirically describes a range 
of dark energy models (Weller & Albrccht 
2002), and permits a first disproof of the 
cosmological constant scenario, if w\ 7^ 
and/or w 7^ — 1. We note that the alterna- 
tive parameterization w{z) = w + w a (l — a) 7 
where a = (1 + z)^ 1 is the usual cosmological 
scale factor, encodes a more physically rea- 
sonable behaviour at high redshift (Lindcr 
2002). In this paper we compute (wo,w a ) 
constraints for one case. 

3. We place Gaussian priors upon the other 
relevant cosmological parameters, the mat- 
ter density f2 m and the Hubble parameter 
h = H(z = 0)/(100 km s" 1 Mpc" 1 ), rather 
than assuming that their values are known 
precisely. 

2.1. Method Summary 

As described in Paper I, the philosophy of 
our analysis is to maintain maximum indepen- 
dence from models. When measuring the acous- 
tic oscillations from simulated data, we divide out 
the overall shape of the power spectrum via a 
smooth 'reference spectrum'. We then fit a sim- 
ple empirically-motivated decaying sinusoid to the 
remaining oscillatory signal. Hence we do not uti- 
lize any information encoded by the shape of the 
power spectrum. This shape may be subject to 
smooth broad-band systematic tilts induced by 
such effects as complex galaxy biasing schemes, 
quasi-linear growth of structure, a running primor- 
dial spectral index, and redshift-space distortions. 
However, it would be very surprising if any of 
these phenomena introduced oscillatory features 



in fc-space liable to obscure the distinctive acous- 
tic peaks and troughs. We note that any model 
where the probability of a galaxy forming depends 
only on the local density field leads to linear bias 
on large scales (Coles 1993; Scherrer & Weinberg 
1998). Furthermore, linear biassing is observed to 
be a very good approximation on large scales (e.g. 
Peacock & Dodds 1994; Cole et al. 2005). This is 
in agreement with numerical simulations of galaxy 
formation which show that galaxies and/or mas- 
sive haloes faithfully reproduce the acoustic oscil- 
lations (Springel et al. 2005; Angulo et al. 2005). 

Of course, a full power spectrum template 
should be fitted to real data as well: our aim here 
is to derive robust, conservative lower limits to the 
efficacy of baryon oscillations experiments, using 
only the information contained in the oscillations. 

An important point is that the fractional error 
in the measured galaxy power spectrum, crp/P, af- 
ter division by a smooth overall fit, is independent 
of the absolute value of P(k) if the error budget 
is dominated by cosmic variance rather than by 
shot noise. In this sense, an incorrect choice of 
the underlying model power spectrum in our sim- 
ulations docs not seriously affect the results pre- 
sented here. Having secured a detection of the 
acoustic signature, if one is then prepared to model 
the underlying power spectrum - correcting for 
such systematic effects as non-linear gravitational 
collapse, redshift-space distortions and halo bias 
- then more accurate constraints on cosmologi- 
cal parameters would follow (see Eisenstein et al. 
2005). 

In summary (see Section 2.2 for a more detailed 
account): we generate a model matter power spec- 
trum in the linear regime using the fitting formu- 
lae of Eisenstein & Hu (1998), assuming a pri- 
mordial spectral index n — 1 (as suggested by 
inflationary models) and fiducial cosmological pa- 
rameters VL m = 0.3, h = 0.7 and baryon fraction 
^b/^m = 0.15, broadly consistent with the latest 
determinations (e.g. Spergel et al. 2003). In Pa- 
per I we showed that the cosmological constraints 
are fairly insensitive to the exact value of Qb/Q m 
(Figures 7-8 in Paper I; increasing fl^ results in a 
somewhat higher baryonic oscillation amplitude, 
hence a more precise measure of the standard 
ruler). We assume that the shape of P(k) does 
not depend on the dark energy component, and 
take the z = normalization ag = 1. The model 



4 



P(k) is then used to generate Monte Carlo realiza- 
tions of a galaxy survey covering a given geometry, 
deriving redshifts and angular co-ordinates for the 
galaxies using a fiducial flat cosmological constant 
model. The realizations are then analyzed for a 
grid of assumed dark energy models. P(k) is mea- 
sured using a Fast Fourier Transform (FFT) up 
to a maximum value of k determined by a conser- 
vative estimate of the extent of the linear regime 
at the redshift in question (see Paper I, Figure 
1). The measured power spectrum is fitted with a 
decaying sinusoid with the 'wavelength' as a free 
parameter. By comparing the fiducial wavescale 
(determined using the values of £l m h 2 and Q^h 2 
in conjunction with a standard fitting formula, e.g. 
Efstathiou & Bond (1999)) with the distribution 
of fitted wavescales across the Monte Carlo real- 
izations, we can reject each assumed dark energy 
model with a derivable level of significance. We 
assume a flat universe (in Section 2.4, Figure 6 
below we indicate how our dark energy measure- 
ments weaken with declining knowledge of fi^). 

Throughout this paper we ensure that our 
model galaxy surveys contain sufficient objects 
that the contribution of shot noise to the error 
in the power spectrum is sub-dominant to that of 
cosmic variance. In an analytical treatment (e.g. 
Tegmark 1997), the relative contributions of shot 
noise and cosmic variance can be conveniently ex- 
pressed in terms of the quantity n x P, where n 
is the typical number density of galaxies in the 
survey volume and P is the galaxy power spec- 
trum evaluated at some typical scale measured 
by the survey. Analytically, the errors due to 
shot noise and to cosmic variance are equal when 
n x P = 1. For the simulations described in this 
paper, we uniformly populated the survey volume 
with sufficient galaxies that n x P = 3, where 
P is evaluated at a characteristic scale k = 0.2 h 
Mpc -1 . The surface density of galaxies required 
to achieve this condition is illustrated in Figure 1. 

2.2. Detailed fitting methodology 

Here we outline the Monte-Carlo approach we 
have implemented in our code. In fact for most of 
our analysis we utilized a simple approximation as 
explained in Section 2.3: to be specific we omitted 
steps 8 and 9 below, which are very expensive in 
computational resources (we test this approxima- 
tion in the Appendix). 
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Fig. 1. — The surface density of galaxies per unit 
redshift, dN/dz, required to achieve power spec- 
trum shot noise levels nxP = 1 (dashed line; equal 
error contribution due to shot noise and to cosmic 
variance) and n x P = 3 (solid line; adopted in our 
simulations). The power spectrum P is evaluated 
at a 'typical' scale k = 0.2 h Mpc -1 ; at zero red- 
shift this value is taken to be P = 2500 h~ 3 Mpc 3 
(and is scaled to higher redshifts using the square 
of the linear growth factor D\(z)). If the galax- 
ies possess a linear bias b, then the amplitude of 
their power spectrum scales as b 2 and a propor- 
tionately lower surface density is required. For 
6 = 1 galaxies in a redshift slice of width Az = 0.2 
at z = 1, a density ss 1800 deg -2 is required to 
achieve n x P = 3. 
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1. A fiducial cosmology is chosen for the sim- 
ulation: for example (O m , h, wo, wi) = 
(0.3,0.7,-1,0). A fiducial baryon fraction 
is selected (fib/^m = 0.15). 

2. A survey redshift range (z m ; n , z max ) and 
solid angle is specified: for example 1.0 < 
z < 1.3 and 1000 deg 2 . The survey sky ge- 
ometry is assumed to be bounded by lines 
of constant right ascension and declina- 
tion of equal angular lengths. The three- 
dimensional geometry is therefore 'conical' 
and the convolving effect of this window 
function is included. 

3. Using the fiducial cosmology, a cuboid for 
FFTs is constructed whose sides (L x ,L y , L z ) 
are just sufficient to bound the survey vol- 
ume. Note, only the enclosed cone is popu- 
lated by galaxies in order that the window 
function effect is treated properly. 

4. A model matter power spectrum P maS s(k, z — 
0) is computed for the chosen parame- 
ters (Q m ,Qb,h) from the fitting formula 
(Eiscnstein & Hu 1998), assuming a z = 
normalization (78 = 1 and a primordial 
power-law slope n = 1. The survey slice 
is assumed to have an 'effective' redshift 
z e tt = (zmin + z max )/2. The power spectrum 
is scaled to this redshift using the linear 
growth factor Di(z), obtained by solving 
the full second-order differential equation 
(see e.g. Linder & Jenkins 2003) to enable 
us to treat non-ACDM cosmologies: 

Pgal(Meff) = Pm,ss(k,0)D 1 (z cS ) 2 b 2 (1) 

where we use a constant linear bias factor b 
for the clustering of galaxies with respect to 
matter. The value b = 1 is assumed for our 
surveys, unless otherwise stated. 

5. A set of Monte Carlo realizations (number- 
ing 400 for all simulations presented here) is 
then performed to generate many different 
galaxy distributions consistent with P ga i(k), 
as described in steps 6 and 7. 

6. A cuboid of Fourier coefficients is con- 
structed with grid lines set by dki = 2^/ Li, 
with a Gaussian distribution of amplitudes 



determined from P ga i(fc), and with random- 
ized phases. The gridding is sufficiently fine 
that the Nyquist frequencies in all directions 
are significantly greater than the smallest 
scale for which a power spectrum is ex- 
tracted (i.e. the linear/non-linear transition 
scale at the redshift z c s). 

7. The Fourier cuboid is FFTed to determine 
the density field in the real-space box. This 
density field is then Poisson sampled within 
the survey 'cone' to determine the number 
of galaxies in each grid cell. 

8. Using the fiducial cosmological parameters, 
this distribution is converted into a simu- 
lated catalogue of galaxies with redshifts and 
angular positions, for each Monte Carlo re- 
alization. 

9. We now assume a trial cosmology: for exam- 
ple (Q m ,/i,u;o,wi) = (0.3,0.7,-0.9,0). The 
co-moving co-ordinates of the galaxies are 
computed in the trial cosmology as would 
be done by an observer without knowledge 
of the true cosmology. 

10. The power spectrum of the simulated sur- 
vey for the trial cosmology is measured using 
standard estimation tools (Feldman, Kaiser, 
& Peacock 1994). Power spectrum modes 
in Fourier space are divided into bins of 
(fcperp, fcpar) where, if the cc-axis is the radial 
direction, fc par = k x and fc 2 orp = k 2 + k 2 z . 

11. An error bar is assigned to each power spec- 
trum bin using the variance measured over 
the Monte Carlo realizations. Note that the 
distribution of realizations also encodes any 
covariances between different power spec- 
trum bins, although the scale of correla- 
tions in fe-space is expected to be very small 
(compared to the separation of the acous- 
tic peaks) for the very large survey volumes 
considered here. 

12. The measured P(fc per p, fcpar) is divided by a 
smooth 'reference spectrum' following Paper 
I (i.e. the 'no-wiggles' spectrum of Eisen- 
stein & Hu 1998), and the result is fitted 
with a simple empirical formula, modified 
from Paper I to permit separate fitting of the 
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sinusoidal scale in the radial and tangential 
directions: 




where k 2 — fc pcrp + k 2 aT . The free parame- 
ters are then the tangential and radial sinu- 
soidal wavescales (Ap Cr p, A par ) together with 
the overall amplitude A. 

We can now assign a probability to the trial 
cosmology. The Monte Carlo realizations produce 
a distribution of fitted wavescales A per p and A par 
(Figure 2) for the trial cosmology. Using these 
trial cosmological parameters, we can determine 
the length of the characteristic ruler Athcory using 
a standard fitting formula for the sound horizon 
integral (Equation 1, Paper 1) in terms of f2 m , 
n h and h (Efstathiou & Bond 1999). (We re- 
mind the reader that Athcory is set in the early 
Universe and is insensitive to dark energy param- 
eters.) The location of the value of Athcory in the 
distribution of tangential (radial) wavescales over 
the Monte Carlo realizations allows us to assign 
a probability for the trial cosmological parame- 
ters. For example: if Athcory lies at the 16 th per- 
centile of the distribution, the rejection probability 
is 2 x (50 - 16) = 68%. Note that the simulated 
observer does not need to know the fiducial cos- 
mological parameters (including the dark energy 
model) to perform this analysis with real data. 

2.3. The streamlined approach 

The full methodology outlined above is too 
computationally intensive for the exploration 
of a full grid of trial cosmological parameters 
(Cl m ,h, Wo,Wi). In practice we pursue a stream- 
lined approach that adopts some simple approx- 
imations. In the Appendix we use a test case to 
demonstrate that the results are equivalent to the 
utilization of the full methodology 

In our streamlined approach, we exploit the fact 
that the accuracy of measurement of A pcrp is a very 
good approximation of the precision with which we 
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Fig. 2. — Distribution of fitted tangential and ra- 
dial wavescales for 400 Monte Carlo realizations 
of a survey covering volume lOVsioan and redshift 
range 0.75 < z < 1.25. The solid circle marks the 
expected value of the acoustic scale and the dashed 
lines bound regions containing 68% of the fitted 
values. The tangential wavescale may be deter- 
mined more accurately than the radial wavescale 
because more tangential modes are available (i.e. 

a cylindrical shell of radius fc porp = ^Jk 2 + k 2 cen- 
tred upon the radial (x-) axis). 
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can recover the quantity x(z e ff)/s, where x{z e s) 
is the co-moving radial distance to the effective 
redshift of the survey and s is the value of the 
sound horizon at recombination. This is because 
(in the flat-sky approximation) the value of x con- 
trols physical tangential scales in the slice (as dis- 
placements Ar are governed by Ar = x x Aff). 
Similarly, the accuracy of measurement of A par 
is equivalent to that of x'(z e ff)/s, where x'(z) = 
dx/dz — c/H(z) (since Ar = x' x Az). The value 
of s appears in the denominators because a sys- 
tematic shift in the standard ruler scale implies a 
similar variation in the recovered physical scales x 
and dx/dz: cosmic distances are measured in units 
of the sound horizon at recombination (equiva- 
lently, we may think of this measuring rod as the 
distance to the CMB: s = Dqmb x 6a, where 9a 
is the angular scale separating the CMB acoustic 
peaks). 

Therefore, rather than re-constructing the 
galaxy distribution using a trial cosmology, we 
instead fitted wavelengths A per p and A par directly 
in the fiducial cosmology (i.e. omitting steps 8 
and 9 above). The 68% scatter in these fits across 
the Monte Carlo realizations was assigned as a 
la Gaussian error in the values of x(z e s)/s and 
x'(z e g)/s, respectively. The likelihood contours 
for the trial cosmological parameters were then 
deduced using standard expressions for dx/dz and 
x in terms of (fi m , h, w n ,wi) and a fitting formula 
for s in terms of (f2 m , f2b, h) (see Section 2.4). 

This streamlined approach assumes that: 

1. The scatter in fitted wavelengths is indepen- 
dent of the values of the cosmological pa- 
rameters. In detail, changing the cosmo- 
logical parameters will alter the cosmic vol- 
ume surveyed between z m i D and z max , and 
therefore the errors in the recovered power 
spectrum in any bin (and hence the accu- 
racy with which the sinusoidal scale may be 
determined). However, these variations can 
be neglected for small perturbations about a 
fiducial cosmology. 

2. The values of [x(z c g) 7 x'(z e g)] control tan- 
gential and radial scales, respectively. This 
statement is exact for a flat sky, and holds 
approximately for the conical geometry as- 
sumed here (if the survey solid angle is not 
too large). 



These approximations are tested in the Ap- 
pendix. 

2.4. Likelihoods for dark energy models 

The procedure thus far has permitted us to re- 
cover values and statistical errors for the quanti- 
ties x(z c s)/s and x'(z c s)/s for each survey red- 
shift bin. These measurements are statistically 
independent to a good approximation (this is evi- 
denced by the distribution of fitted wavelengths in 
Figure 2 being close to an ellipse aligned parallel 
with the axes, displaying only a weak tilt). Fig- 
ure 3 illustrates the simulated recovery of x(z)/s 
and x'(z)/s for some Monte Carlo realizations of 
a 1000 deg 2 survey, in redshift slices of width 
Az = 0.5 from z = 0.5 to z = 3.5. Results for 
the accuracy of recovery of x(z)/s and x'(z)/s in 
these redshift slices are listed in Table 1. 

We parameterize the dark energy model using 
an equation-of-state w(z) = wq + w\z. The ac- 
curacies of x(z)/s and x'(z)/s are then used to 
infer joint constraints over a grid of (wo,Wi), us- 
ing the standard formulae for x'(z) = c/H(z) and 
x(z) = J* x'(z')dz', where the Hubble constant 
H{z) at redshift z is a function of (f2 m , h, Wo, w\). 
For the sound horizon s we used the formulae 
of Efstathiou & Bond (1999, equations 18-20) in 
terms of (f2 m , fib, h) so we effectively assume that 
the effect of dark energy at early times is insignif- 
icant. We explore a range of uncertainties for 
these quantities below. For each grid point in the 
(wq , u>i ) plane we derive a likelihood for each red- 
shift slice by marginalizing over Gaussian priors 
for f2 m and Q m h 2 (the natural variables - see be- 
low) centered upon f2 m = 0.3 and h = 0.7. The 
sound horizon is only weakly dependent on fib, 
therefore we simply fixed the value Vl^h 2 = 0.022 
(we checked that the likelihood contours remained 
unchanged for reasonable variations in tl^h 2 ). The 
overall (68%, 95%) likelihood contours were then 
determined by multiplying together the individ- 
ual likelihoods inferred from the measurements of 
x(z)/s and x'(z)/s for each redshift slice. 

Figure 4 displays the resulting (wo,w\) con- 
tours for a 10,000 deg 2 survey. A further ap- 
proximation has been used to generate this plot 
(also tested in the Appendix): that the accura- 
cies of the fitted wavelengths determined for the 
1000 deg 2 simulation (Figure 3) may be scaled 
by a factor \/T0. This is simply equivalent to 
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Fig. 3. — Simulated measurements of x(z)/s and x'(z)/s in six redshift slices for the first four Monte Carlo 
realizations of a 1000 deg 2 survey. The plotted values in each redshift bin are inferred from the fractional 
deviation of the fitted tangential and radial wavelengths from their fiducial values and the error bars represent 
the standard deviation across the realizations. The various model curves represent the fractional deviations 
of x(z) and x'{z) as the dark energy model is varied from the fiducial point (vjq, w±) = (—1,0). The plots 
illustrate that these data would be sufficient to rule out dark energy parameters (—0.8, 0) and (—1, 0.4) with 
high confidence; but, owing to the cosmic degeneracy between wq and wi, a model (wq,Wi) — (—0.8,-0.4) 
cannot be confidently excluded. Note that the evolution of w(z) — Wq + w\z is cut off at z = z cut = 2. 
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Survey 


Area 


z- 


-bin 


Accuracy 


Accuracy 




(deg 2 ) 






x/s (%) 


x'/s (%) 


spec- z 


1000 


0.5 


- 1.0 


2.7 


3.8 






1.0 


- 1.5 


1.4 


2.1 






1.5 


- 2.0 


1.2 


2.0 






2.0 


- 2.5 


1.0 


1.9 






2.5 


-3.0 


1.2 


2.0 






3.0 


-3.5 


1.1 


1.9 


spec- z 


10000 


0.5 


- 0.7 


1.7 


2.7 






0.7 


-0.9 


1.1 


2.0 






0.9 


- 1.1 


1.0 


1.5 






1.1 


- 1.3 


0.7 


1.4 






1.3 


- 1.5 


0.6 


1.2 


KAOS 


1000 


0.5 


- 1.3 


1.6 


2.6 




400 


2.5 


- 3.5 


1.2 


2.3 


photo- z do = 0.03 


2000 


0.5 


- 1.5 


2.3 








1.5 


- 2.5 


1.4 








2.5 


-3.5 


1.3 





Table 1: Simulated precision of recovery of the quantities x(z e s)/s and x'(z e s)/s from a series of future 
spectroscopic and photometric redshift surveys. For spectroscopic redshift surveys, these quoted accuracies 
(for area Ai) may be approximately scaled to other survey areas (A 2 ) by multiplying by a factor \f~A\fA/i 
(since the errors in the power spectrum measurement 5P oc 1/y/V). In order to scale the photometric 
redshift results to surveys with different r.m.s. redshift scatters 02, we can multiply by a further factor 
\fa~2 1 txi (since the number of usable Fourier modes m scales as 1/cr, and SP oc 1/^/m). These simple 
scalings will break down (1) if the redshift range changes significantly, owing to the changing position of the 
non-linear transition, and (2) in the regime where we are just resolving the oscillations, when improvement 
is better than y/V. The measurement precision of x/s and x'/s is determined by Monte Carlo realizations 
and is accurate to ± 0.1%. 
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sub-dividing the 10,000 deg 2 survey into 10 sep- 
arate independently-analyzed pieces. If we do 
not make this additional approximation, unfeasi- 
bly large Fourier transforms are required to handle 
the size of the survey cuboid. Furthermore, the 
statistical independence of x(z e ff) and x'(z e g) in a 
given redshift slice becomes weaker, as this inde- 
pendence rests upon the flat-sky approximation. 
We note that: 



1 . There is a significant degeneracy in each red- 
shift slice between wq and wi, because ap- 
proximately the same cosmology is produced 
if wq becomes more negative and wi becomes 
more positive. The axis of degeneracy is 
a slow function of redshift, which improves 
this situation somewhat as we combine dif- 
ferent redshift slices. 

2. As redshift increases, the radial oscillations 
provide decreasingly powerful constraints 
upon the dark energy model, because H(z) 
becomes independent of (wq ,wi). This con- 
clusion is valid if we are perturbing about 
the cosmological constant model (—1,0), but 
will not be true in general for models with 
wi =/= for which dark energy may affect 
dynamics at higher redshift. 

3. The tightness of the likelihood contours 
in the (wq,wi) plane depends significantly 
upon the fiducial dark energy model (see 
Figure 14). As wq and loi become more 
positive, dark energy grows more influential 
at higher redshifts and the simulated sur- 
veys constrain the dark energy parameters 
more accurately, despite the fact that the 
surveyed cosmic volume is decreasing. 

When generating Figure 4, we assume that the 
values of f2 m (= 0.3) and h (= 0.7) are known per- 
fectly. Although there is some useful cancellation 
between the trends of the distance scale x(z) and 
the standard ruler scale s with the value of O m /i 2 , 
there is nevertheless some residual dependence of 
the experimental performance on the accuracy of 
our knowledge of both £l m and h, independently. 

In this paper we choose not to combine our re- 
sults with cosmological priors from specific pro- 
posed experiments (as could be achieved by com- 
bining Fisher matrix information, for example). 



0.5 < z < 1.0 ( 

1.0 < z < 1.5 ( 

1.5 < z < 2.0 ( 

2.0 < z < 2.5 ( 

2.5 < z < 3.0 ( 

3.0 < z < 3.5 ( 

Combined (1ct,2ct) 



0.5 < z < 
10000 deg' 
Aw Q = 0.027 
Aw, = 0.060 




-1.1 



Fig. 4. — Likelihood contours for dark energy 
model parameters (wo,wi) for a 10,000 deg 2 sur- 
vey. Contours are shown for six individual redshift 
slices ranging from z — 0.5 to z — 3.5 (68%), to- 
gether with the combined result (68%, 95%). For 
this plot, precise knowledge of Sl m (= 0.3), h 
(= 0.7) and f\ (= 0) is assumed. In Figures 5 
and 6 we plot sets of contours for different priors 
in fl m , h and fik- The quoted errors Awq and Ami 
are produced by marginalizing over the other dark 
energy parameter (i.e. wi and Wo, respectively). 
This plot is obtained by scaling the inferred mea- 
surement accuracies of x/s and x /s from a 1000 
deg 2 simulation by a factor \/l0. 
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We prefer to keep our presentation in general 
terms by marginalizing over different priors for Q nl 
and h. We recognize that this approach will not 
capture all the parameter degeneracies inherent 
in future CMB, large-scale structure or supernova 
surveys, but nevertheless we can robustly quan- 
tify the accuracy of knowledge required for these 
other cosmological parameters such that their un- 
certainties do not dominate the resulting error 
in dark energy parameters. Combinations with 
any specific future experiment can be achieved by 
using our results listed in Table 1 which repre- 
sent the fundamental observables recovered by this 
method: x{z) and x'{z) in units of the sound hori- 
zon. 

In general terms, one degree of freedom in other 
parameters is constrained by the excellent mea- 
surement of the physical matter density fl m h 2 af- 
forded by the CMB angular power spectrum: ac- 
curacies of about 3% {a(Q m h 2 ) w 0.004) and 1% 
(a{tt m h 2 ) fa 0.001) are possible with the WMAP 
and Planck satellites, respectively (Balbi et al. 
2003). However, a second independent constraint 
on a combination of f2 m and h is also required. 

This is illustrated by Figure 5, in which we 
re-compute the overall likelihoods in the (wq,w-l) 
plane, marginalizing over the WMAP and Planck 
errors in fl m h 2 together with a second indepen- 
dent Gaussian prior on fi m . We conclude that for 
a survey of 10,000 deg 2 , £l m must be known with 
an accuracy a(Q m ) ~ 0.01 (in conjunction with 
the WMAP or Planck determination of fl m h 2 ) in 
order that this uncertainty is not limiting. Note 
that as the cosmic volume surveyed increases, the 
prior requirements of knowledge of Q m become 
more stringent. For a 1000 deg 2 experiment, only 
c(O m ) < 0.03 is required (see Figure 10). 

Large-scale structure and/or supernovae con- 
straints together with the first-year WMAP CMB 
data currently deliver cr(f2 m ) ~ 0.02-0.04 as- 
suming a flat Universe (e.g. Spergel et al. 2003; 
Tegmark et al. 2004; Cole et al. 2005). Balbi et 
al. (2003) quote a(h) ~ 0.02 as attainable with 
Planck (equivalent to cr(O m ) ~ 0.02), even al- 
lowing for uncertainty in the dark energy model. 
However, in light of Figure 5, this may not be suf- 
ficient for a 10,000 deg 2 baryon oscillations survey. 

Table 2 lists some 68% confidence ranges for 
dark energy parameters (wq , w\ ) for a range of 
survey configurations and cosmological priors, as- 
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Fig. 5. — Likelihood contours (68%) for dark 
energy model parameters (wq,wi) for the same 
(combined) survey as Figure 4, considering vari- 
ous different Gaussian priors upon the values of 
tt m h 2 and fl m , with standard deviations a(fl m h 2 ) 
and a(fl m ). We assume fit = 0. For the values 
of a(fl m h 2 ) we make assumptions representative 
of the WMAP and Planck CMB experiments, for 
each of these cases we further consider cr(f2 m ) = 
0.01 and 0.03. For a survey of 10,000 deg 2 cover- 
ing redshift range 0.5 < z < 3.5, WMAP suffices 
to determine the value of Q m h 2 , but we require 
an independent constraint cr(f2 m ) ~ 0.01 in order 
that this uncertainty is not limiting. 
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Survey 


Configuration 


a(r» m /i 2 ) 


<7(Q m ) 


<t(w ) 


o-(tui) 


spec- z 


(10000 deg 2 , 0.5 < z < 3.5) 








0.03 


0.06 






WMAP 


0.01 


0.04 


0.07 






WMAP 


0.03 


0.08 


0.10 






Planck 


0.03 


0.07 


0.09 


KAOS 


(1000 deg 2 , 0.5 < z < 1.3) + (400 deg 2 , 2.5 < z < 3.5) 








0.17 


0.48 






WMAP 


0.03 


0.27 


0.63 






WMAP 


0.05 


0.34 


0.71 




{z ~ 1) + (1000 deg 2 , 1.5 < z < 2.5) 








0.10 


0.26 


SKA 


(20000 deg 2 , 0.5 < z < 1.5) 








0.04 


0.11 






Planck 


0.01 


0.05 


0.13 






Planck 


0.03 


0.11 


0.18 


photo-z 


(10000 deg 2 , 0.5 < z < 3.5, a = 0.01) 








0.07 


0.19 






WMAP 


0.01 


0.20 


0.57 






WMAP 


0.03 


0.30 


0.95 




(10000 deg 2 , 0.5 < z < 3.5, er = 0.03) 








0.12 


0.32 






WMAP 


0.01 


0.23 


0.66 






WMAP 


0.03 


0.31 


0.95 




(2000 deg 2 , 0.5 < z < 3.5, cr = 0.01) 








0.19 


0.51 






WMAP 


0.01 


0.25 


0.72 






WMAP 


0.03 


0.31 


0.95 



Table 2: Simulated 68% confidence ranges for the dark energy parameters (wo,wi) for a series of future 
spectroscopic and photometric rcdshift surveys, assuming a fiducial cosmology (—1, 0). A range of (Gaussian) 
priors on the values of Q m h 2 and Q m are considered. The WMAP and Planck measurement precisions of 
a(Q m h 2 ) are assumed to be 0.004 and 0.001, respectively. We assume Ok = 0. 
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suming a fiducial model (—1, 0). 

Throughout this paper we assume a spatially- 
flat (S7k = 0) cosmology. In Figure 6 we compute 
how the likelihood contours in the (wq,wi) plane 
for our 10,000 deg 2 survey degrade as our knowl- 
edge of fik weakens. We note that current deter- 
minations of the curvature (<r(Ok) ^ 0.02, Spergel 
et al. 2003; Eisenstein et al. 2005) are almost ad- 
equate for this proposed experiment. Of interest 
is Bernstein (2005) which noted that the the com- 
bination of baryon oscillations with weak lcnsing 
constraints leads to direct breaking of degenera- 
cies of curvature with dark energy and allows Ok 
to be measured without any assumptions about 
the equation of state. 
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Fig. 6. — Likelihood contours (68%) for dark 
energy model parameters (wq,wi) for the same 
(combined) survey as Figure 4, considering var- 
ious different Gaussian priors upon the value of 
Ok = with standard deviations cr(Ok). We con- 
sider both holding fixed the values of O m and h, 
and marginalizing over these parameters. 

2.5. Comparison with the Fisher matrix 
methodology 

It is worth comparing our methodology and re- 
sults with the Fisher matrix approach utilized by 
Seo & Eisenstein (2003) for the simulation of bary- 
onic oscillations experiments. The input data and 
assumptions arc not entirely consistent in the two 
cases, but we can make a reasonably direct com- 
parison of, for example, Figure 5 in Seo & Eisen- 
stein (2003) with Figure 9 in this study. In such 
comparisons we find that the accuracies of deter- 
mination of (wq, Wi) are consistent within a factor 



of about 1.5, with the Fisher matrix method yield- 
ing tighter contours. 

This is not surprising: the Fisher matrix 
method uses information from the whole power 
spectrum shape (which will also be distorted in an 
assumed cosmology, e.g. by the Alcock-Paczynski 
effect) whereas in our approach, this shape is di- 
vided out and only the oscillatory information 
is retained. The bulk of the Fisher information 
does appear to originate from the sinusoidal fea- 
tures (see Figure 5 in Hu & Haiman (2003) and 
Section 4.4 of Seo & Eisenstein (2003)). How- 
ever, the improvement in dark energy precision 
resulting from fitting a full power spectrum tem- 
plate to the data may be up to 50% compared 
to our 'model-independent' treatment. Our re- 
sults should provide a robust lower limit on the 
accuracy, as intended. 

We also note that the Fisher approach provides 
the minimum possible errors for an unbiased es- 
timate of a given parameter based upon the cur- 
vature of the likelihood surface near the fiducial 
model. As a result the projected error contours for 
any combination of two parameters always form an 
ellipse (e.g. Figure 5 of Seo & Eisenstein (2003)). 
In our approach we explore the whole parameter 
space and estimate probabilities via Monte Carlo 
techniques: the error contours are thus larger and 
not necessarily elliptical (e.g. our Figure 9). 

A further difference between the appearances 
of our Figure 9 and Figure 5 in Seo & Eisen- 
stein (2003) is a noticeable change in the tilt of 
the principal degeneracy direction in the (wq, w±) 
plane. The reason for this is readily identified: Seo 
& Eisenstein additionally incorporate the CMB 
measurement of the angular diameter distance to 
recombination into their confidence plots. We 
choose not to do this in order to expose the 
low redshift independent cosmological constraints 
from galaxy surveys and to isolate our results 
from the effects of any unknown behaviour of the 
equation of state (which could extend beyond our 
wo,wi formalism) between z ~ 4 and z ~ 1100. 

In summary, we are encouraged by the rough 
agreement in values of Awq and Awi between 
these two very different techniques. They rep- 
resent respectively more conservative/robust and 
best possible dark energy measurements from fu- 
ture surveys for baryon oscillations. 
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3. Detectability and accuracy of wavescale 
extraction 

In this Section we take a step back from ques- 
tions of dark energy and employ our simulation 
tools to re-consider the fundamental question of 
the detectability of the acoustic oscillations in 
P(k) as a function of survey size and redshift cov- 
erage. 

The oscillations in the galaxy power spectrum 
are a fundamental test of the paradigm of the ori- 
gin of galaxies in the fluctuations observed in the 
early Universe via the CMB. 

Detection of the oscillations would be an ex- 
tremely important validation of the paradigm. 
Furthermore our standard ruler technique cannot 
be confidently employed unless the sinusoidal sig- 
nature in the power spectrum can be observed 
with a reasonable level of significance. 

Recently, analysis of the clustering pattern of 
SDSS Luminous Red Galaxies at z = 0.35 (a vol- 
ume of ~ 3.5Vsioan) has yielded the first convinc- 
ing detection of the acoustic signal and application 
of the standard ruler (Eisenstein et al. 2005). Al- 
though this survey does not have sufficient redshift 
reach to strongly constrain dark energy models, 
this result is an important validation of the tech- 
nique. Analysis of the final database of the 2dF 
Galaxy Redshift Survey has also yielded some vi- 
sual evidence for baryonic oscillations (Cole et al. 
2005). Here we make the distinction between de- 
tection of oscillations and detection of a baryonic 
signal =/= in the clustering pattern: baryons 
produce an overall shape distortion in P(k) as well 
as the characteristic pattern of oscillations. In this 
Section we define the 'wiggles detectability' as the 
confidence of rejection of a smooth 'no wiggles' 
model (i.e. the best-fitting smooth reference spec- 
trum of step 12 in Section 2.2). This is a different 
quantity to the '3.4-sigma confidence' of observ- 
ing Sib 7^ reported by Eisenstein et al. (2005). 
However, the techniques roughly agree on the mea- 
surement accuracy of the standard ruler (see the 
discussion of SDSS LRGs in our Paper I, Figure 
3). 

Figure 7 tracks the detection significance 
against the percentage accuracy of recovery of the 
standard ruler for surveys covering three different 
redshift ranges: 0.25 < z < 0.75, 0.75 < z < 1.25 
and 2.75 < z < 3.25. The lines connect points sep- 



arated by survey volume intervals of lT^sioan; the 
solid circles denote volumes (2, 5, 10) Vsioan- F°r 
the purposes of this Section, an angle-averaged 
(isotropic) power spectrum is used rather than a 
power spectrum separated into tangential and ra- 
dial components. We also only consider the pure 
vacuum ACDM model: the detectability is primar- 
ily driven by the cosmic volume surveyed, which is 
a relatively slow function of dark energy parame- 
ters. The detection significance is calculated from 
the average over the Monte Carlo realizations of 
the relative probability P Ic i of the smooth 'no- 
wiggles' model and best-fitting 'wiggles' model, 
where 

P rc l = CXp [-(Xno-wig ~ Xbcst-wig)/ 2 ] (3) 



O 



o 

e o 

0) ,- 
u 

CL - 

'Z> d 

D 

-Q _ 

B o 

°- d 



: 1 1 1 1 1 i : 


- 1a _ 


r 2a * - y 












r / ' f J* 












i / / 












3(7 *' / 












'•1 L 

■ 1 ] 
r ! / * 












; 4u j\J 1 


0.25 


< 


z 


< 


0.75 


1 t 1 

r : ' 


- - 0.75 


< 


z 


< 


1.25 1 














1 1 


2.75 


< 


z 


< 


3.25 


: 5a : 
■ ■ ■ i ■ 



-13.5 1 2 5 10 



Accuracy 6A/X (per cent) 

Fig. 7. — Variation of the 'wiggles detectability' 
(i.e. the significance of rejection of a no-wiggles 
model), and the accuracy with which the char- 
acteristic sinusoidal scale may be extracted, with 
survey volume. Three surveys are considered at 
different redshifts, centred upon z — 0.5, z = 1 
and z = 3. The lines join points separated by vol- 
ume intervals of lVsi oan ; the three solid circles on 
each line denote volumes (2, 5, 10) Vsioan- 

We note that: 

1. In order to obtain a significant (3a) rejec- 
tion of a no- wiggles model without using any 
power spectrum shape information, several 
Vsioan must be surveyed. For the surveys 
centred at z = (0.5, 1, 3) the required volume 
(in units of Vsioan) is approximately (7, 5, 5). 
For the redshift ranges listed above, this 
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corresponds to survey areas (2300, 700, 500) 
deg 2 . 

2. For a fixed wiggles detection significance, the 
accuracy of recovery of the standard ruler 
increases with redshift. This is due to the 
larger available baseline in k at higher red- 
shift, owing to the extended linear regime. 
Two full oscillations are visible at z = 1; 
four are unveiled by a survey at z = 3. 

3. For a fixed wavescale accuracy, the detection 
significance decreases with redshift, because 
the amplitude of the oscillations is damped 
with increasing k. As noted above, at higher 
redshifts there are more acoustic peaks avail- 
able, thus a less significant measurement of 
each individual peak may be tolerated. 

4. The distribution of values of P ro i (equation 
3) over the Monte Carlo realizations is sig- 
nificantly skewed (see the discussion in Blake 
& Bridle (2005)). The median value of P rol 
represents a more confident detection than 
the average plotted in Figure 7. 

Figure 8 displays some Monte Carlo power spec- 
trum realizations of three surveys: (0.5 < z < 1.3, 
1000 deg 2 ), (2.5 < z < 3.5, 400 deg 2 ) and (0.5 < 
z < 1.5, 10,000 deg 2 ). The total volumes mapped 
in units of Vsioan are (10, 8, 133), respectively. The 
total numbers of galaxies observed in each survey 
(to ensure nxP ~ 3) are (6, 2, 85) x 10 6 (assuming 
linear bias b = 3 for the z <~ 3 survey and 6=1 
otherwise). The first two redshift surveys could 
be performed by a next-generation wide-field op- 
tical spectrograph such as the KAOS instrument 
proposed for the Gemini telescopes (Bardcn ct al. 
2004, http://www.noao.edu/kaos/). The third 
survey is possible in 6 months using the Square 
Kilometre Array to detect HI emission line galax- 
ies (Abdalla & Rawlings 2004) or from a space 
mission (see Section 4.4). 

4. Dark energy measurements from realis- 
tic spectroscopic redshift surveys 

We now consider the prospects of performing 
these experiments with realistic galaxy redshift 
surveys. As noted above and in Paper I, such sur- 
veys must cover a minimum of several hundred 
deg 2 at high redshift, cataloguing at least several 



hundred thousand galaxies, in order to obtain sig- 
nificant constraints upon the dark energy model. 

4.1. Existing surveys 

These requirements are orders of magnitude 
greater than what has been achieved to date. 
Some existing surveys of high-redshift galaxies are 
the Canada- France Redshift Survey (CFRS; a few 
hundred galaxies covering ~ 0.1 deg 2 to z m 1.3; 
Lilly et al. 1995), the survey of z ~ 3 Lyman- 
break galaxies by Steidel et al. (2003) (roughly 
a thousand galaxies across a total area w 0.4 
deg 2 ). Most other high redshift spectroscopic sur- 
veys (e.g. Abraham et al. 2004; Cimatti et al. 2002; 
Steidel et al. 2004) cover equally small areas < 1 
deg 2 . 

Some larger surveys are in progress: the 
DEEP2 project (Davis et al. 2003), using the 
DEIMOS spectrograph on the Keck telescope, 
aims to obtain spectra for 60,000 galaxies (3.5 
deg 2 , 0.7 < z < 1.4); the VIRMOS redshift survey 
(Le Fevre et al. 2003), using the VIMOS spec- 
trograph at the VLT, will map 150,000 redshifts 
over 16 deg 2 (considering the largest-area compo- 
nent of each). Neither of these existing projects 
comes close to meeting our goals, primarily due to 
the limitations of existing instrumentation. The 
spectrographs used to perform these surveys have 
typical fields-of-view (FOV) of diameter 10 - 20' 
and are unable to cover hundreds of deg 2 in a 
reasonable survey duration. 

4.2. New ground-based approaches (opti- 
cal/IR) 

Some proposed new optical instrumentation ad- 
dresses this difficulty, permitting spectroscopic 
exposures over considerably larger FOVs using 
the 8-metre telescopes that arc required to ob- 
tain spectra of sufficient quality at these red- 
shift depths. For example, the KAOS project 
for the Gemini telescopes (Barden et al. 2004, 
http://www.noao.edu/kaos/) is a proposal for 
a 1.5 deg diameter FOV, 4000 fibre- fed optical 
spectrograph. There are two proposed redshift 
surveys: 900,000 (0.5 < z < 1.3) galaxies over 
1000 deg 2 , and 600,000 (2.5 < z < 3.5) galax- 
ies across 400 deg 2 . These surveys would together 
take ~ 170 clear nights using the 8-metre Gemini 
telescope with realistic exposure times computed 
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Fig. 8. — Power spectrum realizations of three example high-redshift surveys, varying survey volume and 
redshift range. In all cases, we plot the angle-averaged power spectrum divided by the smooth reference 
spectrum. The dashed curve is the theoretical input P(k) and the solid line is the best fit of our simple 
decaying sinusoidal function (Paper I, equation 3). The x-axis is marked in units of k (in h Mpc^ 1 ) and 
represents the extent of the linear regime at the redshift in question (i.e. z c g = (z m i n + z max )/2). The rows 
of the Figure represent surveys with geometries (0.5 < z < 1.3, 1000 deg 2 ), (2.5 < z < 3.5, 400 deg 2 ) and 
(0.5 < z < 1.5, 10,000 deg 2 ); the columns display the first three Monte Carlo realizations in each case. 
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for the KAOS instrument sensitivity. The redshift 
ranges are driven by the strong spectral features 
available for redshift measurement in optical wave- 
bands in a relatively short exposure time. The 
z ~ 1 range is cut off at z = 1.3 by the [Oil] 
emission line and the calcium H & K lines shifting 
to red/infrared wavelengths > 0.9 fim where the 
airglow is severe and conventional CCD detectors 
have low efficiency; the z ~ 3 component is driven 
by observing Lya in the blue part of the optical 
range. 

The w(z) measurements resulting from the 
proposed KAOS surveys, computed using the 
methodology of Section 2, are displayed in Fig- 
ure 9 (see also Tables 1 and 2). We show both 
the z ~ 1 and z ~ 3 contributions separately, and 
the joint constraint. We assume linear bias factors 
b = (1,3) for the z ~ (1,3) simulations, respec- 
tively. The measurement precision of the dark 
energy parameters is Awo ~ 0.2 and Awi » 0.4, 
significantly better than current supernovae con- 
straints (Riess et al. 2004). In statistical terms 
the KAOS performance is somewhat poorer than 
that projected for the proposed SNAP mission 
(Aldering et al. 2002), but the acoustic oscillations 
method is significantly less sensitive to errors of a 
systematic nature. Figure 9 illustrates that mod- 
els with w\ < are harder to exclude owing to 
the diminishing sensitivity of cosmic distances to 
dark energy in this region of parameter space. 

We note that the KAOS measurement of H (z « 
3) from the radial component of the z ~ 3 power 
spectrum provides little information about the 
dark energy model if we are perturbing around 
the cosmological constant: at z ~ 3, the dynam- 
ics of the Universe are entirely governed by the 
value of fl m h 2 , and H(z) is almost independent 
of (wo,Wi). However, the value of x(z rs 3) in- 
ferred from the tangential component of the z ~ 3 
power spectrum does depend on dark energy, be- 
cause x(z = 3) is an integral of dx/dz = c/H(z) 
from z = to z — 3, which is influenced by dark 
energy at lower redshifts. The z ~ 3 constraint 
thus reduces to a significant degeneracy between 
wq and itfi, as observed in Figure 9, although mod- 
els with w± > can still be ruled out by this red- 
shift component. The degeneracy is less severe for 
the z ~ 1 component owing to the availability of 
both H(z) and x{z) information. 

Figure 9 assumes that we have perfect prior 
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Fig. 9. — Measurement of a dark energy model 
w(z) = wq + w\z using a simulated survey with 
the KAOS spectrograph consisting of two compo- 
nents, 0.5 < z < 1.3 (1000 deg 2 ) and 2.5 < z < 3.5 
(400 deg 2 ). Constraints are displayed for each 
redshift component (68%), and for both surveys 
combined (68%, 95%). Note that the likelihoods 
are generated over a much wider (wq,wi) space 
than displayed in the Figure. The solid circle de- 
notes the fiducial cosmology, (w ,Wi) — (—1,0). 
There is a significant degeneracy between wq and 
W\, particularly for the z ~ 3 survey if w± < 0, ow- 
ing to the lack of sensitivity of H(z — 3) to dark 
energy in this case. For this plot, perfect knowl- 
edge of fi m (= 0.3), h (= 0.7) and ilk (= 0) is as- 
sumed. The standard deviations quoted for Awo 
and Awi (i.e. half the interval between the 16 th 
and the 84 th percentiles) result from marginalizing 
over the other parameter (i.e. w\ and Wq, respec- 
tively). 
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knowledge of £l m (= 0.3) and h (= 0.7). Figure 
10 relaxes this assumption, illustrating the effect 
of including Gaussian priors upon fi m and fl m h 2 
of various widths. 
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Fig. 10. — Likelihoods contours (68%) of dark en- 
ergy parameters (wo, wi) for the same KAOS red- 
shift surveys as Figure 9, marginalizing over the 
cosmological parameters (fi m , f2 m /i 2 ). We assume 
a WMAP Gaussian prior upon fl m h 2 (as in Fig- 
ure 5) and consider a range of different Gaussian 
priors for f2 m . In each case, standard deviations 
are quoted for wq and w\ as before. We assume 
n k = 0. 

A drawback of this survey design is the absence 
of the redshift range 1.5 < z < 2.5, sometimes 
called the 'redshift desert'. There are no strong 
emission lines accessible to optical spectrographs 
in this interval: existing surveys of this region have 
required very long exposure times to secure spec- 
tra (Abraham et al. 2004), but this could be reme- 
died by near-IR or near-UV spectroscopy of bright, 
star-forming galaxies (Steidel et al. 2004). We now 
consider the usefulness of these additional obser- 
vations as regards measuring dark energy, along 
with the observational practicalities. 

First, we investigate the effect of this red- 
shift range upon measurements of the dark en- 
ergy model w(z) = wq + w\z (assuming precise 
prior knowledge of Q m and h for the purposes of 
this comparison; Figure 10 indicates how accu- 
rately these parameters must be known in order 
that their uncertainty is not limiting). In Fig- 
ure 11, we remove the z = 3 component of the 
proposed KAOS experiment and extend the lower- 
redshift 1000 deg 2 survey across the redshift range 



1.5 < z < 2.5, divided into two independent slices 
of width Az = 0.5. The likelihood constraints 
in (u>o, w\) space tighten appreciably, by a further 
factor of two, principally due to the 1.5 < z < 2.0 
component, for which H{z) still yields useful infor- 
mation about (wq, w\). In Figure 12 we add back 
in the z = 3 data; the dark energy measurements 
do not significantly improve. We conclude that 
coverage of the redshift desert would be highly de- 
sirable if it could be achieved. Simply increasing 
the area of the z ~ 3 component does not help 
nearly as much (Awi is improved by ~ 25% for 
1000 deg 2 at z ~ 3). This simplistic analysis is of 
course no substitute for a proper survey optimiza- 
tion assuming fixed total time or cost. 
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Fig. 11. — Dark energy measurements result- 
ing from 1000 deg 2 surveys spanning the redshift 
ranges (0.5 < z < 1.3) and (1.5 < z < 2.5). We 
assume perfect prior knowledge of (O m , h, fik)- 

We note, however, that there are other persua- 
sive scientific reasons to include a z — 3 survey 
component (Eisenstein 2002), amongst them: 

1. A galaxy redshift survey at z = 3 unveils 
the linear power spectrum down to unprece- 
dentedly small scales k w 0.5 h Mpc -1 , mea- 
suring linear structure modes that cannot be 
accessed using the CMB. 

2. If dark energy is insignificant at z = 3, then 
measurement of the acoustic oscillations in 
such a redshift slice enables the standard 
ruler to be calibrated in a manner indepen- 
dent of the CMB. 

3. Our present analysis assumes that the fidu- 
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Fig. 12. — Dark energy measurements result- 
ing from 1000 deg 2 surveys spanning the redshift 
ranges (0.5 < z < 1.3) and (1.5 < z < 2.5), to- 
gether with a 400 deg 2 component covering (2.5 < 
z < 3.5). We assume perfect prior knowledge of 
(tt m ,h, O k ). 

cial dark energy model is a cosmological con- 
stant, (wq,w±) = (—1,0). Based on current 
data, we have very little information about 
the value of Wi (for the latest supernova 
analysis of Riess et al. (2004), u(wi) « 0.9). 
Should u>i 7^ 0, the influence of w(z) upon 
higher-redshift dynamics could become more 
significant. A general redshift survey opti- 
mization, which is beyond the scope of this 
paper, should address the range of w{z) to 
be explored (Bassett et al. 2005). 

Figure 13 considers a different parameterization 
for dark energy w(z) = wq + w a (l — a) = wq + 
w a /(l + z) (Linder 2002), where a = (1 + z)" 1 is 
the usual cosmological scale factor, which encapsu- 
lates a more physically realistic behaviour at high 
redshifts z > 1 in comparison to w(z) = wq + w\z. 
The rate of change of w with redshift in the two 
models is dw/dz — w\ = w a /(l + z) 2 < w a hence 
for a given survey we expect the size of the error 
for w a to exceed that for w%. Figure 13 illustrates 
this for the same survey configuration as Figure 
12. 

Figure 14 uses the KAOS surveys plus a 1.5 < 
z < 2.0 extension to illustrate how the tightness of 
the likelihood contours in the (wq , w\ ) plane is a 
strong function of the fiducial dark energy param- 
eters, as discussed in Section 2.3. We computed 
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Fig. 13. — Dark energy measurements resulting 
from the same surveys as Figure 12, using the al- 
ternative parameterization w{z) = wq + w a (l — a). 
We assume perfect prior knowledge of (f2 m , h, Q^). 
Note that the likelihoods are generated over a 
much wider (wo,w a ) space than displayed in the 
Figure. 

the linear growth factor for non-ACDM models 
by solving the appropriate second-order differen- 
tial equation (e.g. Linder & Jenkins 2003). If 
wi > then dark energy is more significant at 
higher redshifts and the model parameters can 
be constrained more tightly, despite both the re- 
sulting decrease in the available cosmic volume 
in a given redshift range and the movement of 
the linear/non- linear transition to larger scales 
(smaller k). Note that our cut-off to the evolu- 
tion of w(z) — wo + w\z at z cut = 2 ensures that 
dark energy does not dominate at high redshift for 
our model with w\ > 0. 

Next, we consider the practicalities of observing 
galaxies in the redshift range 1.5 < z < 2.5 (the so- 
called 'optical redshift desert'). Considering near 
infra-red wavebands first: the Ha 6563 A line is 
accessible in the l-2^m band over the interval 
0.5 < z < 2. This regime is the non-thermal infra- 
red, in which the sky is sufficiently dark to permit 
high-redshift spectroscopy (e.g. the emission-line 
observations of Glazebrook et al. (1999) and Pet- 
tini et al. (1998)). The critical opportunity offered 
here is that the Ha emission line is potentially very 
bright at redshifts z > 1, owing to the steep evo- 
lution in the star-formation rate of galaxies in the 
Universe over < z < 1 (Hopkins, Connolly, & 
Szalay 2000). 
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Fig. 14.— The tightness of the (68%, 95%) likeli- 
hood contours in the (wo,Wi) plane is a strong 
function of the fiducial cosmological model. In 
this Figure we illustrate, for three different fidu- 
cial models (wq 7 wi) = (—1,0), (—0.8,0.4) and 
(—1.2,-0.4), dark energy measurements result- 
ing from 1000 deg 2 surveys spanning the redshift 
ranges (0.5 < z < 1.3) and (1.5 < z < 2.0), to- 
gether with a 400 deg 2 component covering (2.5 < 
z < 3.5). We assume perfect prior knowledge of 
{Sl m ,h, Ok)- 

Let us consider a redshift slice 1.5 < z < 
1.7. Our requirement for shot noise sampling 
(n x P = 3) translates into a required surface 
density (3900/6 2 ) deg -2 within this slab (Figure 
1), where b is the linear bias factor of the sur- 
veyed galaxies. The luminosity function of Ha 
emitters at z > 1 has been reasonably well deter- 
mined from NICMOS slitless grism surveys using 
the Hubble Space Telescope (Hopkins, Connolly, 
& Szalay 2000; Yan et al. 1999). Using the Hop- 
kins et al. luminosity function, we find that we 
need to reach a line flux limit of 1.1 x 10 -16 ergs 
cm -2 s _1 in order to reach the required surface 
density (assuming 6=1). The z = 1 Ha lumi- 
nosity function has apparently evolved strongly 
in comparison with z — (Gallego et al. 1995) 
but the measurements are fairly robust: we can 
double-check the luminosity function determina- 
tions by simply counting objects in the NICMOS 
surveys above our required line flux limit. In the 
Hopkins, Connolly & Szalay sample there are 4 
galaxies in the redshift range 1.5 < z < 1.7 above 
this flux limit, spread over 4.4 arcmin 2 , yielding a 
surface density ~ 3300 deg -2 . The Ha identifica- 



tions are also very reliable: the Yan et al. sample 
was observed in optical wavebands by Hicks et al. 
(2002), confirming > 75% of the Ha identifications 
via associated [Oil] emission at the same redshift. 
This agrees with expectations: analytic models of 
evolving line emission show that Ha should domi- 
nate at these flux levels at 1-2/im over other lines. 

This is encouraging, because these bright lines 
are accessible in relatively modest exposures. Let 
us assume an 8-metre telescope, a 25% efficient 
R = 4000 near-IR spectrograph and detector 
(with a dispersion of 5 A per arcsec), and consider 
an object observed in an aperture of size 0.8 arcsec 
x 0.8 arcsec, covering 2x2 pixels and containing 
half the light (the Yan et al. objects have com- 
pact half-light radii of 0.2-0.7 arcsec). The signal- 
to-noise ratio is determined by detector readout 
noise, sky background and dark current. We will 
assume a readout noise of 4 electrons (this is the 
stipulated requirement of the James Webb Space 
Telescope ( JWST) detectors) and an inter-OH sky 
background of 1.2 photons s -1 ran -1 arcsec" 2 m~ 2 
in the H-band (the OH night sky line forest is well- 
resolved at R = 4000) which is a typical value 
measured at Gemini observatory 1 . We will ne- 
glect dark current, which is equivalent to assum- 
ing that this dark current is much lower than the 
sky counts. Given these assumptions, our Ha flux 
limit at z = 1.6 corresponds to a signal-to- noise 
ratio of 20 in a 600 second integration, and we 
find the observation is background limited for any 
readout noise < 15 electrons. 

This exposure time is encouragingly short, and 
implies that using a 1 degree FOV spectrograph, 
one could survey 1000 deg 2 in only 20 nights. Ad- 
mittedly our instrument specification is optimistic, 
especially for readout noise; however, it can be re- 
laxed considerably whilst still achieving exposure 
times < 1 hour. Assuming a fibre spectrograph 
we estimate that observing 3900 objects simulta- 
neously at this resolution would only require 2-3 
detector arrays of size 4096 x 4096 to cover the H- 
band. Obviously more objects or broader wave- 
length coverage would require more detectors or 
time. Finally we note that, at least in princi- 
ple, the exposure times are sufficiently short that 
one could imagine performing such a survey on a 

See http:/ /www. gemini.edu/sciops/ObsProcess/obsConstraints/ 
ocSkyBackground.html 
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smaller-aperture (4-metre class) telescope. 

The potential problem with the approach out- 
lined above is that it is not known 'a priori' which 
of the galaxies identified in deep images will be 
Ha bright, or the redshifts of these galaxies. Pos- 
sibly these data could be successfully predicted 
from other information (e.g. broad-band colours or 
sub- mm/radio fluxes), but this may be unreliable. 
Targetting fainter Ha fluxes resulting from realis- 
tic target selections would obviously require longer 
exposure times. Further work on this problem is 
required in order to determine the distinguishing 
properties of known Ha-bright galaxies. 

A second potential difficulty is the effect of the 
night-sky OH emission lines in potentially mak- 
ing inaccessible certain redshift ranges, in a com- 
plex pattern. Since these redshift ranges are very 
narrow, the result is the removal of a series of 
redshift spikes at known locations in the radial 
window function. In order to assess the likely 
consequences, we manufactured a synthetic n(z) 
possessing narrow gaps where no galaxies could 
be observed, i.e. when the redshifted Ha emis- 
sion line coincided with an OH line or landed in 
the water-absorption hole between the J and H 
bands. We assumed a spectrograph operating at 
a resolution R = 4000, which implied that 68% of 
the 1.1 < z < 1.7 redshift interval was accessible. 
With our simulation tools we then recovered P(k) 
using this n(z) (employing an FFT with sufficient 
gridding in the radial direction to resolve these 
narrow spikes) . The principal effect of the window 
function is to damp the amplitude of the acoustic 
oscillations slightly in the radial direction, leav- 
ing the tangential modes largely unaffected. The 
fractional error (Afc^/Zc^) with which the acous- 
tic scale is recovered is increased by no more than 
25%, which is mostly due to the smaller effective 
survey volume owing to the absence of many thin 
redshift shells. We conclude that the OH lines 
are not a factor which will significantly hamper 
ground-based surveys of the acoustic peaks. 

Bright star-forming galaxies can alternatively 
be observed by targetting the [Oil] 3727A emis- 
sion line in optical wavebands. High-resistivity 
CCD detectors under development can maintain 
quantum efficiency out to 1/xm (Stover et al. 1998) 
corresponding to z = 1.7. The typical Ha: [Oil] in- 
tensity is 2:1; high spectral resolution is again re- 
quired to observe between the OH night sky lines, 



in which case the exposure times would be short 
if one could pre-select the [OII]-bright population. 
This could be more problematic than with Ha as 
there is considerable extra scatter introduced in to 
the line ratio [Oil] /Ha due to variations in metal- 
licity and extinction (Jansen, Franx & Fabricant 
2001). We also note that [Oil] is accessible in the 
non-thermal IR up to z = 5 and therefore is a po- 
tential probe of high-redshift acoustic oscillations. 

The high star-formation rate at earlier cos- 
mic epochs also implies a considerable boost in 
the rest-frame UV luminosity of galaxies. High- 
altitude sites such as Mauna Kea have an atmo- 
spheric cutoff further into the near-UV, which can 
be exploited by blue-optimized spectrographs. For 
example, Steidel et al. (2004) have obtained spec- 
tra of star-forming galaxies over 1.4 < z < 2.5 
using exposure times of only a few hours at the 10- 
metre Keck telescope, reaching a lowest observed- 
frame wavelength of 3200A; this corresponds to 
Lya at z = 1.6, or CIV at z = 1.1. The spectral 
region between Lya and CIV is rich in interstel- 
lar lines and ripe for determination of accurate 
redshifts. It is possible that a UV approach may 
be superior to a ncar-IR approach targetting Ha: 
the observed surface density of Steidel et al. sam- 
ple is sufficiently high for our requirements, but 
we note that a UV-optimized design would prob- 
ably require a wide-field slit spectrograph because 
conventional fibres considerably attenuate the UV 
light for long runs (> 20m). 

4.3. New ground-based approaches (ra- 
dio) 

Next-generation radio interferometer arrays, 
such as the proposed Square Kilometre Array 
(SKA; http://www.skatelescope.org; planned 
to commence operation in about 2015), will have 
sufficient sensitivity to detect the HI (21cm) tran- 
sition of neutral hydrogen at cosmological dis- 
tances that are almost entirely inaccessible to cur- 
rent radio instrumentation. This will provide a 
very powerful means of performing a large-scale 
redshift survey: once an HI emission galaxy has 
been located on the sky, the observed wavelength 
of the emission line automatically provides an ac- 
curate redshift. 

The key advantage offered by a radio telescope 
is that it may be designed with an instantaneous 
FOV exceeding 100 deg 2 (at 1.4 GHz), vastly sur- 
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passing the possibilities of optical spectrographs. 
Equipped with a bandwidth of many 100 MHz, 
such an instrument could map out the cosmic web 
(probed by neutral hydrogen) at an astonishing- 
rate: the SKA, if designed with a large enough 
FOV, could locate ~ 10 9 HI galaxies to redshift 
z ss 1.5 over the whole visible sky in a timescale 
of — 1 year (Abdalla & Rawlings 2004; Blake et 
al. 2004). Deeper pointings could probe the HI 
distribution to z ~ 3 over smaller solid angles. 
A caveat is that the HI mass function of galaxies 
has been determined locally (Zwaan et al. 2003) 
but is currently very poorly constrained at high 
redshift. However, for a range of reasonable mod- 
els, the number densities required to render shot 
noise negligible may be attained for HI mass limits 
larger than the break in the mass function (Ab- 
dalla & Rawlings 2004). Another requirement of 
interferometer design is that a significant fraction 
of the collecting area must reside in a core of di- 
ameter a few km, to deliver the necessary surface 
brightness sensitivity for extended 21cm sources. 
Sharp angular resolution is not a pre-requisite, as- 
suming that the observed galaxies are not con- 
fused. 

Figure 15 displays measurements of the dark 
energy model resulting from a 20,000 deg 2 neutral 
hydrogen survey over the redshift range 0.5 < z < 
1.5, analyzing acoustic peaks in redshift slices of 
width Az = 0.2. We note that a smaller 21cm 
survey could be performed in the nearer future 
by SKA prototypes such as the HYFAR proposal 
(Bunton et al. 2003), which may cover several 
thousand deg 2 over a narrower bandwidth (cor- 
responding to 0.8 < z < 1.2). 

4.4. Space-based approaches 

An interesting alternative approach is to use a 
space-based dispersive but slitless survey to pick 
out emission-line objects directly over a broad red- 
shift range. In space, the 1-2/im background is 
1000 times lower (in a broad band) than that ob- 
served from the ground (based upon the JWST 
mission background simulator at distance 3 au 
from the Sun; Petro, Kriss & Stockman 2002); in 
the background-limited regime this gain is equiva- 
lent to a 1000-fold increase in collecting area. The 
dispersing element could be either a large objec- 
tive prism or a grism. The NICMOS surveys men- 
tioned above already demonstrate that this tech- 
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Fig. 15. — Likelihoods (68%) of dark energy pa- 
rameters (wo,u>i) for a future 21cm radio survey, 
using the Square Kilometre Array to map HI emis- 
sion galaxies over 20,000 deg 2 , covering a redshift 
range 0.5 < z < 1.5. Such a survey is possible in a 
timescale of less than 1 year if the SKA is designed 
with a sufficiently large FOV (~ 100 deg 2 at 1.4 
GHz) and bandwidth (~ 200 MHz). We marginal- 
ize over the cosmological parameters (f2 m , f2 m /i 2 ) 
using a range of different Gaussian priors. Our as- 
sumed prior for fl m h 2 is representative of that ob- 
tained by the Planck satellite. We assume f^k = 0. 

nique is possible, but they lack somewhat in FOV 
and spectral resolution compared to what is desir- 
able. 

As an illustrative example, let us consider a 
0.5-metre space telescope with a 60% overall sys- 
tem efficiency working in slitless dispersed imag- 
ing mode, again targeting 0.5" galaxies at redshift 
z = 1.6. We will assume the JWST- like back- 
ground and that the slitless spectra are de-limited 
by a 2000A-wide blocking filter. Just considering 
the sky background, the signal-to-noise ratio in a 
1800s exposure is 5 for our canonical flux limit of 
1.1 x 10~ 16 ergs cm -2 s _1 . This signal-to-noise 
ratio is independent of spectral resolution for un- 
resolved lines much brighter than the continuum 
(the spectral resolution must of course be suffi- 
cient for determining accurate redshifts). Highly- 
dispersive IR materials such as silicon could en- 
able an objective prism approach with R — 500, 
using a simple prime focus imaging system. Such 
a satellite, if equipped with a 1 deg FOV and sen- 
sitive over 1-2/im, could perform a spectroscopic 
survey over the redshift interval 0.5 < z < 2 (using 
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Ha), covering an area of 10,000 deg 2 in a 4 year 
mission, obtaining dark energy constraints very 
similar to those presented in Figure 4. Access- 
ing higher rcdshifts would be possible by either 
extending the wavelength range or going fainter 
with the [Oil] line, in both cases requiring a larger 
diameter mirror. Ambiguous line identification is 
a potential problem, this could be remedied by 
cross-matching with a photometric-redshift imag- 
ing survey. We refer the reader to Glazebrook et 
al. (2004) for a more detailed discussion of such a 
dedicated 'Baryonic Oscillation Probe'. 

5. Dark energy measurements from realis- 
tic photometric redshift surveys 

We next explore the potential of surveys based 
on photometric redshifts for detecting the acoustic 
oscillations and placing constraints upon any cos- 
mic evolution of the dark energy equation-of-state 
w(z). For a detailed treatment of the significance 
of 'wiggles detection' and the accuracy of measure- 
ment of the standard ruler with photometric red- 
shift surveys, we refer the reader to Blake & Bridle 
(2005). Here we present a summary and discuss 
the consequences for dark energy measurement. 

A photometric redshift ('photo-z') is obtained 
from multi-colour photometry of a galaxy: the ob- 
ject is imaged in several broad-band filters, rang- 
ing from the UV to the near-IR, producing a rough 
spectral energy distribution (SED). This observed 
SED is then fitted by model galaxy SEDs as a 
function of redshift to construct a likelihood dis- 
tribution with redshift; the peak of the likelihood 
function indicates the best-fitting redshift. A clas- 
sic example of this approach is the analysis of 
the Hubble Deep Field North (Fernandez-Soto, 
Lanzetta, & Yahil 1999). 

There are obviously many options concerning 
the number and width of filter bands, and their 
placement in the UV-NIR range. Generally at 
least five broad bands are used, and IR coverage is 
essential for constraining galaxies with 1.2 < z < 
2.2 (Bolzonella, Miralles, & Pello 2000). Some ap- 
proaches have used as many as 17 broad + narrow- 
band filters (e.g. Wolf et al. 2003). Many differ- 
ent techniques have been proposed for deriving 
the photo-z (e.g. Csabai et al. 2000; Lc Borgne & 
Rocca-Volmerange 2002; Collister & Lahav 2004). 
The important question for our study is: what is 



the accuracy of the photo-z estimates? 

These errors can be divided into two main 
types. First, there is the random statistical er- 
ror due to noise in the flux estimates and to the 
coarseness of the SED. Typically this is specified 
by a parameter a where 

z 

<7n = t, 7 ~ constant (4) 

(1 + z) w 

and o~ z is the standard deviation of the redshift 
z. We note that cto is approximately constant be- 
cause it is proportional to the spectral resolution 
A/ A A of the set of filters. Chen et al. (2003) ob- 
tained CT = 0.08; the COMB017 survey achieved 
(Jo = 0.03. In a theoretical study, Budavari et al. 
(2001) demonstrated that an optimized filter set 
produced results within the range cto = 0.02 — 0.05, 
depending on the shape of the SED. In general, 
redder galaxies deliver more accurate photo-z's be- 
cause the model colours change faster with red- 
shift. 

The second source of photo-z error is the pos- 
sibility of getting the redshift grossly wrong, ei- 
ther because the set of colours permit more than 
one redshift solution, or because the model SEDs 
are not sufficiently representative of real galaxies. 
Different authors disagree about the magnitude 
of this effect, which depends on the specific filter 
sets, photometric accuracy, spectroscopic calibra- 
tion and photo-z methods used. A useful theoret- 
ical discussion is given in Bolzonella, Miralles, & 
Pello (2000). Blake & Bridle (2005) analyze vari- 
ous 'realistic' redshift error distributions including 
outliers and systematic offsets. 

For our purposes we will ignore systematic er- 
rors and parameterize photo-z performance using 
the value of er alone. A realistic survey will con- 
tain additional systematic redshift errors, thus we 
will obtain lower limits. We will also assume that 
er is a constant, whereas for a realistic survey it 
will depend somewhat on redshift and galaxy type. 

What is the effect of a statistical redshift er- 
ror ((T ) on the measured power spectrum? This 
is fairly easy to estimate analytically. As dis- 
cussed in Paper I (see also Seo & Eisenstein 2003, 
Section 4.5) this photo-z error represents a ra- 
dial smearing of galaxy positions. For example, 
(Jo = 0.03 for a z = 1 galaxy corresponds to an 
error a x w 100 ft -1 Mpc in the radial co-moving 
coordinate. Since smoothing (i.e. convolution) 
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by a Gaussian function in real space is equiva- 
lent to multiplication by a Gaussian function in 
Fourier space, we can model the photo-z effect as 
a multiplicative damping of the 3D power spec- 
trum P(k x , k y , k z ) by a term cxp (— fc 2 cr 2 ) (where 
we choose the x-axis as the radial direction). In 
the following, we implement a flat-sky approxima- 
tion and presume that there is no tangential effect 
(i.e. parallel to the (y, z)-plane). 

Prior to the smearing effect of photometric red- 
shifts, the available Fourier structure modes in the 
linear regime comprise a sphere in Fourier space 
of radius k < fcimear s=s 0.2 h Mpc -1 (where the 
value of linear depends on redshift). Afterwards, 
the damping term exp (— fc 2 cr 2 ) implies that only a 
thin slice of this sphere with \k x \ < 2/a x w 0.02 h 
Mpc -1 is able to contribute useful power spectrum 
signal. This is illustrated schematically in Fig- 
ure 16. Considering those modes contributing to 
a Fourier bin centred about scale k, the reduction 
in usable Fourier space volume corresponds to a 
factor (k x a x /2). At a scale k = 0.2 /i Mpc -1 , 
this represents a loss by a factor of 10 (cr = 0.03, 
2=1). For the same surveyed area, we can hence 
expect the error ranges in the derived power spec- 
trum to worsen by a factor w \/T0 in comparison 
to a spectroscopic survey (although also note that 
the scaling of the error with k also changes from 
<5P sp cc oc k~ x to cLPphoto oc k^ 1 / 2 ). The errors on 
the resulting dark energy parameters will be worse 
by a factor of ~ \/20 than in a similar sized spec- 
troscopic survey once one also accounts for the loss 
on the radial dimension. 

Critically, the radial damping due to photomet- 
ric redshifts results in the loss of any ability to 
detect the acoustic oscillations in the radial com- 
ponent of the power spectrum (in the above exam- 
ple, the power spectrum is significantly suppressed 
for modes with \k x \ > 2/a x w 0.02 ft, Mpc -1 , i.e. 
the whole regime containing the oscillations). We 
are only able to apply the standard ruler in the 
tangential direction. As noted earlier, the radial 
component provides a direct measure ofH(z), con- 
tributing significantly to the measurements of the 
dark energy model. Thus in the above example 
((To = 0.03, z = 1) the resulting errors in the dark 
energy parameters will worsen by a factor closer 
to "\Z20- Alternatively, one could compensate by 
increasing the survey area (and hence the density 
of states in Fourier space) by the corresponding 



factor. We note that with sufficient filter coverage 
and for special classes of galaxy, the photometric 
redshift error er ma y also be reduced to improve 
dark energy performance. 

We simulated photo-z surveys using analogous 
Monte Carlo techniques to those described in Sec- 
tion 2 (sec Blake & Bridle (2005) for a more de- 
tailed account). We introduced a radial Gaussian 
smearing 

<y x = CT (1 + Zeff) x'(z e s) (5) 

into our Poisson-sampled density fields (for a sur- 
vey slice at redshift z e s). When measuring the 
power spectrum we restrict ourselves to modes 
with \k x \ < 2/a x (where the factor of 2 was de- 
termined empirically to be roughly optimal). The 
residual damping in the shape of P(k) is divided 
out using the known Gaussian damping expres- 
sion. We bin the power spectrum modes in accor- 
dance with the total length of the Fourier vector 
k = ^Jk x + k y + fc|, noting that only tangential 

modes (with k x w 0) are being counted. We fit a 
ID decaying sinusoid (Paper I, equation 3) to the 
result. The scatter of the fitted wavescales across 
the Monte Carlo realizations is interpreted as the 
accuracy of measurement of the quantity x(z e s)/s 
(given that only tangential modes are involved). 
We uniformly populated our survey volumes such 
that nx P = 3. 

We consider two different photo-z accuracies: 
(To = 0.03, representing the typical fidelity of the 
current best photo-z studies, and Co = 0.01, which 
we somewhat arbitrarily adopt as an upper limit 
to future improvements. Blake & Bridle (2005) 
consider a much wider range of possibilities. We 
note that in our methodology, reducing the value 
of (To by some factor is equivalent to covering a 
proportionately larger survey area. Figure 17 dis- 
plays some Monte Carlo realizations of measured 
power spectra for these photometric redshift sur- 
veys, assuming a redshift range 0.5 < z < 1.5. We 
assume a survey solid angle of 10,000 deg 2 , but 
also consider a smaller project (2000 deg 2 with 
(T = 0.01). 

Figure 18 illustrates the resulting measurement 
of the dark energy parameters (wq,Wi) for these 
survey configurations, assuming that we can span 
the redshift range 0.5 < z < 3.5 (also see Tables 
1 and 2). These (w ,wi) contours are computed 
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Spectroscopic survey |k|<0.2 



Photo-z survey 
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Fig. 16. — Illustration of the loss of modes in three dimensional Fourier space (k x , k y , k z ) by smearing along 
the x (redshift) axis. Only long-wavelength modes > smearing length are unsuppressed and the |A;| < 0.2 
sphere is truncated to a thin disk. 
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Fig. 17. — Power spectrum realizations for three example photometric redshift surveys, varying survey area 
and photo-z accuracy (parameterized by <7q). In all cases we assume a survey redshift interval 0.5 < z < 1.5 
(i.e. z g = 1). Only radial Fourier modes with \k x \ < 2/a x (where <r x = <Jo (1 + z e g) x'(z e ff)) are binned 
when measuring P(k), which is plotted divided by a smooth reference spectrum. As in Figure 8, the dashed 
curve is the theoretical input P(fc) and the solid line is the best fit of our simple decaying sinusoidal function 
(Paper I, equation 3). The £-axis is marked in units of k (in h Mpc -1 ) and represents the extent of the linear 
regime at redshift z c g. The rows of the Figure represent surveys with parameters (2000 deg 2 , Co = 0.01), 
(10,000 deg 2 , cr = 0.03) and (10,000 deg 2 , <t = 0.01); the columns display the first three Monte Carlo 
realizations in each case. 
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using the same method as Section 2.4, utilizing 
measurements of x(z)/s in three redshift bins of 
width Az = 1 (for real data, narrower redshift 
slices would be used and the results co-added). 
Each of these redshift constraints corresponds to 
a degenerate line in the (wq,Wx) plane (i.e. con- 
strains one degree of freedom in the dark energy 
model) but, as described in Section 2.4, the di- 
rection of degeneracy slowly rotates with redshift: 
the combination of the likelihoods for each redshift 
bin results in closed contours. 
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Fig. 18. — Dark energy measurements resulting 
from three different photometric redshift surveys 
covering the interval 0.5 < z < 3.5. We assume 
perfect prior knowledge of (O m , h, Ok). Table 2 
lists results for a range of different priors. 

Note that for photometric redshift surveys, the 
cosmological priors on fl m h 2 and fl m required to 
achieve a given measurement precision of (wq, w±) 
are much tighter than for spectroscopic surveys 
(see Table 2). This is because in the absence 
of H(z) information, the photo-z survey must 
achieve a significantly tighter measure of x(z) to 
recover a corresponding measurement accuracy of 
the dark energy parameters, rendering it more sus- 
ceptible to uncertainties in f2 m and h. 

The real figure-of-merit for comparison of prac- 
tical instruments is the accuracy with which the 
dark energy model can be measured for a fixed 
total observing time or at a fixed cost. The myr- 
iad details of comparing large imaging cameras 
with large spectroscopic systems are beyond the 
scope of this paper. However, we note that the 
proposed Large Synoptic Survey Telescope (LSST; 
Tyson 2002) could image half of the entire sky to 



the required depth (V ~ 25) in multiple colours 
every 25 nights; such a survey would produce 
dark energy constraints comparable with a spec- 
troscopic (e.g. KAOS) survey of 1000 deg 2 (the 
latter requiring 170 nights) if ctq = 0.01 could be 
achieved. We regard this level of photometric red- 
shift accuracy as unlikely for ground-based sur- 
veys. We note that the KAOS measurements ad- 
ditionally constrain H(z) leading to qualitatively 
more robust measurement of dark energy. Further 
the KAOS survey could be improved by adding 
more area, whereas once the LSST has observed 
the whole celestial hemisphere, there is obviously 
no further gain in w(z) information from further 
passes. However, since 'all-sky' deep multi-colour 
surveys are being performed for other scientific 
reasons (e.g. cosmic shear analysis), it is of con- 
siderable value to utilize these data for baryonic 
oscillation studies. Furthermore, as discussed in 
detail by Blake & Bridle (2005), a photometric 
redshift survey of several thousand square degrees 
may provide constraints competitive with spectro- 
scopic surveys in the short term. 

6. Conclusions 

This study has extended the methodology of 
Blake & Glazebrook (2003) to simulate measure- 
ments of the cosmic evolution of the equation-of- 
state of dark energy w(z) from the baryonic oscil- 
lations, using the simple parameterization w(z) = 
wq + w\z. The methodology used is very similar 
to that of Paper I, treating the primordial bary- 
onic oscillations in the galaxy power spectrum as 
a standard cosmological ruler, whilst dividing out 
the overall shape of the power spectrum in order to 
maximize model-independence. In this study we 
make the improvement of fitting independent ra- 
dial and tangential wavescales, showing that this is 
directly equivalent to measuring Da(z) and H(z) 
in a series of redshift slices in units of the sound 
horizon. This results in improved constraints upon 
{wo,wi). We have tested the approximations en- 
coded in our approach and found them all to be 
satisfactory, increasing our confidence in the in- 
ferred error distributions for the dark energy pa- 
rameters. The simulated accuracies for (wo,Wi) 
are roughly consistent with other estimates in the 
literature based on very different analysis meth- 
ods. 
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Our baseline 'KAOS-like' optical surveys of <~ 
1000 deg 2 , which can be realized by the next gen- 
eration of spectroscopic instruments at ground- 
based observatories, deliver measurements of the 
dark energy parameters with precision Awo f=s 
0.15— 0.2 and Awi w 0.3— 0.4. In statistical terms, 
these constraints are poorer than those which may 
be provided by a future space-based supernova 
project such as the SNAP proposal. We note, 
however, that the baryonic oscillations method ap- 
pears to be substantially free of systematic error, 
with the principal limitation being the amount of 
cosmic volume mapped. In addition, any mea- 
surement of deviations from a cosmological con- 
stant model is of sufficient importance for physics 
that entirely independent experiments would be 
demanded to confirm the new model. A next- 
generation radio telescope with a FOV w 100 deg 2 
at 1.4 GHz, performing a rcdshift survey of 21cm 
emission galaxies over several 1000 deg 2 , may be 
available on a similar timescale. 

We have considered the observational possibil- 
ities of more extensive baryonic oscillation ex- 
periments covering a significant fraction of the 
whole sky. Such a survey (encompassing 0.5 < 
z < 3.5) may be straight-forward using a dedi- 
cated several- year space mission with slitless spec- 
troscopy. In radio wavebands, the Square Kilo- 
metre Array would be able to survey the en- 
tire visible sky out to z w 1.5 in 6 months, if 
equipped with a sufficiently large FOV. These ex- 
periments would deliver extremely precise mea- 
surements of the dark energy model with accu- 
racy Aw w 0.03 - 0.05 and Awi w 0.06 - 0.1 
and would be invaluable to pursue if a significant 
non-vacuum dark energy signal was detected by 
smaller surveys. 

We have also explored in detail the potential 
of photometric rcdshift optical imaging surveys 
for performing baryonic oscillations experiments. 
The loss of the radial oscillatory signal, due to 
the damping caused by the redshift errors, implies 
that we can no longer recover information about 
the Hubble constant at high redshift. However, 
the baryonic oscillations can still be measured us- 
ing tangential Fourier modes. A deep 2000 deg 2 
imaging survey with excellent photometric red- 
shift precision (cr < 0.03) would allow the oscil- 
lations to be detected (2.5<r significance). Useful 
constraints upon the dark energy model are pos- 



sible if ~ 20,000 deg 2 can be surveyed (such an 
experiment with <jq = 0.03 is roughly equivalent 
to a ~ 1000 deg 2 spectroscopic survey). 

We conclude that the baryonic oscillations in 
the clustering power spectrum represent one of the 
rare accurate probes of the cosmological model, 
possessing the potential to delineate cleanly any 
cosmic evolution in the equation-of-state of dark 
energy, via accurate measurements of Da(z) and 
H(z) in a series of redshift slices. Importantly, 
such an experiment is likely to be substantially free 
of systematic error. Recent observations of SDSS 
Luminous Red Galaxies at z = 0.35 have provided 
the first convincing detection of the acoustic sig- 
nature and validation of the technique. The chal- 
lenge now is to create the large-scale surveys at 
higher rcdshifts required for mapping the proper- 
ties of the mysterious dark energy. 
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A. Tests of approximations 

A.l. Effect of redshift-space distortions 

As mentioned in Section 2.1, a power spectrum measured from a real redshift survey is liable to suffer 
systematic shape distortions. We performed a test to verify that we were able to divide out such distortions 
via a smooth fit, prior to fitting for the acoustic oscillations. For a survey covering 0.5 < z < 1.3 and 
1000 deg 2 , we incorporated a redshift-space distortion into each Monte Carlo realization by smearing the 
redshift of each simulated galaxy by the equivalent of a radial velocity of 300 km s _1 . This process has the 
effect of damping the overall power spectrum by a factor increasing with k and amounting to about 20% at 
k = 0.2 h Mpc -1 . After division by the 'no-wiggles' reference spectrum of Eisenstein & Hu (1998), we fitted 
an additional second-order polynomial to the residual, in order to remove the extra shape distortion. We 
found that the scatter of the fitted acoustic scales across the Monte Carlo realizations (i.e. the precision with 
which we could recover the standard ruler) was entirely unchanged by this more complex fitting procedure. 

A. 2. Streamlined vs Full methodology 

In this Section we evaluate the validity of the approximations employed in our 'streamlined method' 
(Section 2.3) with respect to the full methodology (Section 2.2). These approximations are implemented 
because it would otherwise require a prohibitively long computational time to cover the grid of dark energy 
parameters (wo,wi) with sufficient resolution. In particular, we wish to test that: 

1. The constraints upon the dark energy parameters (wo, wi) arising from the fitted radial and tangential 
power spectrum wavescales can be deduced by supposing that these fitted wavescales measure the 
quantities x/s and x 1 /s, respectively. 

2. Dark energy measurements arising from a survey of solid angle A2 may be inferred from those resulting 
from a survey of area A\ < A 2 by multiplying the accuracy of x/s and x'/s deduced from this former 
survey by a factor y/ A1/A2. 

3. A reasonably broad survey redshift interval, z m ; n < z < z max , can be utilized to measure x/s and x'/s 
at an 'effective' redshift z off = (z min + z max )/2. 

In order to test the effect of these approximations, we consider a 'test case' which we analyze using both 
the streamlined and the full method. This test case comprises a spectroscopic survey covering the redshift 
interval 0.7 < z < 1.3 (i.e. approximately the KAOS z <~ 1 component), but with an increased area of 4000 
deg 2 . This is the maximum size we can simulate in a reasonable time. Implementing the full method, we 
first divide the redshift range into three bins of width Az = 0.2. For each bin, we consider a grid of trial 
dark energy parameters (wo, wi) (spaced by Au>o = 0.05, Awi — 0.1) and create 50 Monte Carlo realizations 
for each grid point, deducing the relative probability of each trial cosmology by assessing the position of the 
theoretical standard ruler scale in the distributions of fitted radial and tangential wavescales (as described 
in Section 2.2). We then obtain the overall constraints upon the dark energy model, resulting from the full 
method, by multiplying together these three likelihood maps. We compare these measurements with those 
resulting from the streamlined approach, based upon scaling a single analysis of a (0.7 < z < 1.3, 1000 deg 2 ) 
survey by a factor \/4, and deeming the effective redshift to be z e g = 1. 

These two likelihood maps in the (wq,wi) plane are compared in Figure 19. The grid is necessarily 
quite coarse and noisy for the full methodology, but it is clear that the likelihood patterns are in reasonable 
agreement concerning both width and degeneracy direction. The most noticeable discrepancy appears in the 
region Wq > — 1, W\ < 0, where dark energy has a reduced influence at high redshift but the disagreement is 
of marginal significance. In the regime toi > 0, the agreement is much better. We conclude that within the 
noise of the simulation the two sets of contours are a good match though it would be desirable in the future 
to test this further. 
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Fig. 19. — Comparison of (68%, 95%) likelihood contours deduced using the 'full methodology' (Section 2.2; 
solid lines) and 'streamlined methodology' (Section 2.3; dashed lines) for a test survey covering area 4000 
deg 2 and redshift interval 0.7 < z < 1.3. The two sets of contours are in reasonable agreement. 
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